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OPINION DYNAMICS: INHOMOGENEOUS BOLTZMANN-TYPE 
EQUATIONS MODELLING OPINION LEADERSHIP AND POLITICAL 

SEGREGATION 

BERTRAM DURING^ AND MARIE-THERESE WOLFRAM^ 


Abstract. We propose and investigate different kinetic models for opinion formation, when the 
opinion formation process depends on an additional independent variable, e.g. a leadership or a 
spatial variable. More specifically, we consider: (i) opinion dynamics under the effect of opinion 
leadership, where each individual is characterised not only by its opinion, but also by another 
independent variable which quantifies leadership qualities; (ii) opinion dynamics modelling po¬ 
litical segregation in the ‘The Big Sort’, a phenomenon that US citizens increasingly prefer to 
live in neighbourhoods with politically like-minded individuals. Based on microscopic opin¬ 
ion consensus dynamics such models lead to inhomogeneous Boltzmann-type equations for the 
opinion distribution. We derive macroscopic Fokker-Planck-type equations in a quasi-invariant 
opinion limit and present results of numerical experiments. 


1. Introduction 

The dynamics of opinion formation have been studied with growing attention; in particular in 
the field of physics [151 uni im, in which a new research field termed sociophysics (going back 
to the pioneering work of Galam et al. IE]) emerged. More recently, different kinetic models to 
describe opinion formation have been proposed Ea Ema la mnzi HO]. Such models successfully 
use methods from statistical mechanics to describe the behaviour of a large number of interacting 
individuals in a society. This leads to generalisations of the classical Boltzmann equation for gas 
dynamics. Then the framework from classical kinetic theory for homogeneous gases is adapted to 
the sociological setting by replacing molecules and their velocities by individuals and their opinion. 
Instead of binary collisions, one considers the process of compromise between two individuals. 

The basic models typically assume a homogeneous society. To model additional sociologic 
effects in real societies, e.g. the influence of strong opinion leaders [12], one needs to consider 
inhomogeneous models. One solution, which arises naturally in certain situations (see, e.g. m), 
is to consider the time-evolution of distribution functions of different, interacting species. To some 
extent this can be seen as the analogue to the physical problem of a mixture of gases, where the 
molecules of the different gases exchange momentum during collisions [5]. This leads to systems 
of Boltzmann-like equations for the opinion distribution functions fi = fi{w, t), i = 1,..., n, of n 
interacting species which are of the form 



The Boltzmann-like collision operators Qij describe the change in time of fi{w^t) due to binary 
interaction depending on a balance between the gain and loss of individuals with opinion w. The 
suitably chosen relaxation times Tij allow to control the interaction frequencies. To model exchange 
of individuals (mass) between different species, additional collision operators can be present on 
the right hand side of 0 . which are reminiscent of chemical reactions in the physical situation. 

Another alternative is to study models where the distribution function depends on an additional 
variable as e.g. in [14]. This leads to inhomogeneous Boltzmann-type equations for the distribution 
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function / = f{x,w^t) which are of the following form 

(2) + div:,(^{x,w)f) = ^Q{f,f). 

Clearly, the choice of the field ^ = ^{x,w) which describes the opinion flux plays a crucial role. It 
may not be easy to determine a suitable field from the economic or sociologic problem, in contrast 
to the physical situation where the law of motion yields the right choice. 

In this paper, we give two examples of opinion formation problems, which can be modelled 
using inhomogeneous Boltzmann-type equations. One is concerned with opinion formation where 
the compromise process depends on the interacting individuals’ leadership abilities. The other 
considers the so-called ‘Big Sort phenomena’, the clustering of individuals with similar political 
opinion observed in the USA. Both problems lead to an inhomogeneous Boltzmann-type equation 
of the form Q. 

Our work is based on a homogeneous kinetic model for opinion formation introduced by Toscani 
in [25]. The idea of this kinetic model is to describe the evolution of the distribution of opinion 
by means of mieroseopie interactions among individuals in a society. Opinion is represented 
as a continuous variable w ^ X with X = [—1,1], where ±1 represent extreme opinions. If 
concerning political opinions X can be identified with the left-right political spectrum. Toscani 
bases his model on two main aspects of opinion formation. The first one is a eompromise proeess 
iniiiQiET], in which individuals tend to reach a compromise after exchange of opinions. The 
second one is self-thinking, where individuals change their opinion in a diffusive way, possibly 
influenced by exogenous information sources like the media. Based on both Toscani [25] defines 
a kinetic model in which opinion is exchanged between individuals through pairwise interactions: 
when two individuals with pre-interaction opinion v and w meet, then their post-trade opinions 
u* and w* are given by 

(3) u* = u — ^P{\v — w\){v — w) f}D{v), = w — ^P{\v — w\){w — v) r]D{w). 

Herein, 7 G (0,1/2) is the constant eompromise parameter. The quantities fj and 77 are random 
variables with the same distribution with mean zero and variance a^. They model self-thinking 
which each individual performs in a random diffusion fashion through an exogenous, global access 
to information, e.g. through the press, television or internet. The functions P(-) and P(-) model 
the local relevance of compromise and self-thinking for a given opinion. To ensure that post¬ 
interaction opinions remain in the interval X additional assumptions need to be made on the 
random variables and the functions P(-) and P('), see [25] for details. In this setting, the time- 
evolution of the distribution of opinion among individuals in a simple, homogeneous society is 
governed by a homogeneous Boltzmann-type equation of the form 0 . In a suitable scaling limit, 
a partial differential equation of Fokker-Planck type is derived for the distribution of opinion. 
Similar diffusion equations are also obtained in [23] as a mean field limit the Sznajd model [24] . 
Mathematically, the model in [25] is related to works in the kinetic theory of granular gases [9]. In 
particular, the non-local nature of the compromise process is analogous to the variable coefficient 
of restitution in inelastic collisions [26]. Similar models are used in the modelling of wealth and 
income distributions which show Pareto tails, cf. m and the references therein. 

The paper is organised as follows. We introduce two new, inhomogeneous models for opinion 
formation. In Sectionwe consider opinion formation dynamics which take into account the effect 
of opinion leadership. Each individual is not only characterised by its opinion but also by another 
independent variable, the assertiveness, which quantifies its leadership potential. Starting from 
a microscopic model for the opinion dynamics we arrive at an inhomogeneous Boltzmann-type 
equation for the opinion distribution function. We show that, alternatively similar dynamics can 
be modelled by a multi-dimensional kinetic opinion formation model. In Section [^ we turn to the 
modelling of ‘The Big Sort’, the phenomenon that US citizens increasingly prefer to live among 
others who share their political opinions. We propose a kinetic model of opinion formation which 
takes into account these preferences. Again, the time evolution of the opinion distribution is de¬ 
scribed by an inhomogeneous Boltzmann-type equation. In Section]^ we derive the corresponding 
macroscopic Fokker-Planck-type limit equations for the inhomogeneous Boltzmann-type equations 
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in a quasi-invariant opinion limit. Details on the numerical solvers as well as results of numerical 
experiments are presented in Section Section [^concludes. 

2. Opinion formation and the influence of opinion leadership 

The prevalent literature on opinion formation has focused on election processes, referendums 
or public opinion tendencies. With the exception of 0112 less attention has been paid to the 
important effect that opinion leaders have on the dissemination of new ideas and the diffusion 
of beliefs in a society. Opinion leadership is one of several sociological models trying to explain 
formation of opinions in a society. Certain, typical personal characteristics are supposed to char¬ 
acterise opinion leaders: high confidence, high self-esteem, a strong need to be unique, and the 
ability to withstand criticism. An opinion leader is socially active, highly connected and held in 
high esteem by those accepting his or her opinion. Opinion leaders appears in such different areas 
as political parties and movements, advertisement of commercial products and dissemination of 
new technologies. 

In the opinion formation model in [12] the society is built of two groups, one group of opinion 
leaders, and one of ordinary people, so-called followers. In this model, individuals from the same 
group can influence each others opinions, but opinion leaders are assertive and, although able to 
influence followers, are unmoved by the followers’ opinions. In this model a leader always remains 
a leader, a follower always a follower. Hence it is not possible to describe the emergence or deline 
of leaders. 

In the model proposed in this section we assume that the leadership qualities (like assertiveness, 
self-confidence, ...) of each individual are characterised through an additional independent variable, 
to which we refer short as assertiveness. In the sociological literature the term assertiveness 
describes a person’s tendency to actively defend, pursue and speak out for his or her own values 
preferences and goals. There has been a lot of research on how assertiveness is connected to 
leadership, see for example [ 2 ] and references therein. 

2 . 1 . An inhomogeneous Boltzmann-type equation. Our approach is based on Toscani’s 
model for opinion formation §, but assumes that the compromise process is also influence by 
the assertiveness of the interacting individuals. The assertiveness of an individual is represented 
by the continuous variable x ^ J with J = [—1,1]. The endpoints of the interval i.e. ±1, 
represent a strongly or weakly assertive individual. A leader would correspond to a strongly 
assertive individual. 

The interaction for two individuals with opinion and assertiveness (l’, x) and (re, y) reads as: 

(4) v"" = V — ^C{x^y^v^w){v — w) ^ f}D{v)^ w"" = w — ^C{x^y^v^w){w — v)riD{w). 

The first term on the right hand sides of models the compromise process, the second the 
self-thinking process. The functions C{-) and D{-) model the local relevance of compromise and 
self-thinking for a given opinion, respectively. The constant eompromise parameter 7 G (0,1/2) 
controls the ‘speed’ of attraction of two different opinions. The quantities fj^ y are random variables 
with distribution 0 with variance and zero mean, assuming values on a set S C M. 

Let us discuss the interaction described in @ and its ingredients in more detail. In each such 
interaction, the pre-interaction opinion v increases (gets closer to w) when v < w and decreases 
in the opposite situation; the change of the pre-interaction opinion w happens in a similar way. 
We assume that the compromise process function C can be written in the following form: 

(5) C{x,y,v,w) = R{x,x -y)P{\v -w\). 

In @ we will only allow interactions that guarantee G T. To this end, we assume 

0 < P{\v — rej) <1, 0 < R{x,x — y) <1, 0 < D{v) < 1. 

Let us now discuss useful assumptions on the functions P, R and D. 

As in [25| we assume that the ability to find a compromise is linked to the distance between 
opinions. The higher this distance is, the lower the possibility to find a compromise. Hence, the 
loealisation funetion P(') is assumed to be a decreasing function of its argument. Usual choices 
are P{\v — w\) = l{\v-w\<c} foi* some constant c > 0 and smoothed variants thereof [25] . 
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Figure 1. Choice of function R: in the limit /c ^ oo a step-function is ap¬ 
proached, and the interaction effectively approaches a variant of the leader- 
follower model in [12] (if the assertiveness of individuals were assumed to be 
constant in time). 


On the other hand, we assume the higher the assertiveness level, the lower the tendency of 
an individual to change their opinion. Hence, R{') should be decreasing in its argument, too. A 
possible choice for i?(-) may be 

(6) R{x, x-y) = R{x ~y) = \~\ {k{x - y)) 

for some constant /c > 0. This choice is motivated by the following considerations: Let A and 
B denote two individuals with assertiveness and opinion (t,u) and {y^w)^ respectively. Then the 
particular choice of <11 corresponds to the two cases: 

• If T ~ 1 and ^ ~ — 1, i.e. a highly assertive individual A meets a weakly assertive individual 
B: then R{x — y) ~ 0 (no influence of individual B on A), but R{y — x) ~ 1, i.e. the leader 
A persuades a weakly assertive individual B. 

• If both individuals have a similar assertiveness level and hence \x — y\ is small, there will 
be some exchange of opinion, no matter how large (or small) this assertiveness level is. 

Note that in the limit k ^ oo the choice <© corresponds to an approximation of I — H{z) with 
H denoting the Heaviside function, see Figure In this limit, individuals are either maximally 
assertive in the interaction with a value of the assertiveness variable in ( 0 , 1 ), or minimally assertive 
with assertiveness in ( — 1,0). Effectively, one would recover in the same limit a variant of the 
opinion leader-follower model in m if the assertiveness of individuals were assumed to be constant 
in time. 

We conclude by discussing the choice of D{-). We assume that the ability to change individual 
opinions by self-thinking decreases as one gets closer to the extremal opinions. This reflects the fact 
that extremal opinions are more difficult to change. Therefore, we assume that the self-thinking 
function D{') is a decreasing function of with D{1) = 0 . We also need to choose the set S, i.e. 
we have to specify the range of values the random variables 77 , in can assume. Clearly, it 
depends on the particular choice for D{-). Let us consider the upper bound at re = I first. To 
ensure that individuals’ opinions do not leave X, we need 

u* = u — 7 C(t, 7 /, u, w){v — w) f}D{v) < I 

Obviously, the worst case is re = I, where we have to ensure 

fjD{v) < I — u -h 7 ( 7 ; — I) = (I — 7 ;)(I — 7 ) 

Hence, if D{v)/{I — v) < it suffices to have \fj\ < (I — A similar computation for the 

lower boundary shows that if D{v)/(1 + t;) < K_ it suffices to have \ f}\ < (I — ^)IK_. 

In this setting, we are led to study the evolution of the distribution function as a function 
depending on the assertiveness x G X, opinion w gT and time t G / = /(x, tu, t). In analogy 
with the classical kinetic theory of rarefied gases, we emphasise the role of the different indepen¬ 
dent variables by identifying the velocity with opinion, and the position with assertiveness. In 
this way, we assume at once that the variation of the distribution /{x^wR) with respect to the 
opinion variable w depends on ‘collisions’ between agents, while the time change of distributions 
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with respect to the assertiveness x depends on the transport term. Unlike in physical applications 
where the transport term involves the velocity field of particles, here the transport term contains 
an equivalent ‘opinion-velocity field’ ^ = ^{x^w) which controls the flux depending on the inde¬ 
pendent variables x and w. This is in contrast to the physical situation of rarefied gases where 
the field would simply be given by ^{w) = w. 

The time-evolution of the distribution function / = f{x,w,t) of individuals depending on 
assertiveness x ^ J^ opinion w G X and time t G will obey an inhomogeneous Boltzmann-type 
equation^ 


(7) 


^f{x,w,t) + div^{^{x,w)f{x,w,t)) = tQ[fj)(x,w,t). 


Herein, ^ is the ‘opinion-velocity field’ and r is a suitable relaxation time which allows to control 
the interaction frequency. The Boltzmann-like collision operator Q which describes the change of 
density due to binary interactions is derived by standard methods of kinetic theory. A useful way 
of writing the collision operators is the so-called weak form. It corresponds to consider, for all 
smooth functions 0(re), 

(8) a Q{f, f) (x, w, t)(j){w) dw dx 

= ^ (/ ~ (P{^))f{x,v,t)f{y,w,t)dvdwdx^ , 

where (•) denotes the operation of mean with respect to the random quantities r), 7^. 

Choosing (t){w) = 1 as test function in ^ and denoting the mass by 


p{fm = u /(x, re, t) dw dx 


implies conservation of mass, 

d/?(/)(e 

dt 


= ~j j Q{f,f){x,w,t)dwdx = d. 


If there is point-wise conservation of opinion in each interaction in @ (e.g. choosing (7 = 1 and 

p* -h re* = t; -h re 

or, more generally, conservation of opinion in the mean in each interaction in 0 (e.g. choosing 

C^l), 

(-e* -h re*) = v -\-w, 

then the first moment, the mean opinion, is also conserved. This can be seen by choosing (j){w) = re 
in (§. Denoting the mean opinion by 


um = a /(x, re, t)w dw dx, 


we obtain 


dm(/)(e 

dt 


= iy y Q{f,f){x,w,t)wdwdx = 0. 


We still have to specify the ‘opinion-velocity field’ ^(x,re). A possible choice can be 
(9) ^(x,re) = G{x,w,t){l — x^)^, 

where G = G(x,w,t) models the in- or decrease of the assertiveness level, while the prefactor 
(1 — x^)^ ensures that the assertiveness level stays inside the domain J". 

This leads to an inhomogeneous Boltzmann-type equation of the following form 

7/ + div3;(G(a;,M;,i)(l - A)“/) = 
ot T 


(10) 
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2 . 2 . A multidimensional Boltzmann-type equation. Alternatively, we can consider that 
both, change of opinion and and change of the assertiveness level happen through binary col¬ 
lisions. In this case we have the following interaction rules: 

(11a) = V — ^R{x, X — y)P{\v — w\){v — w) f}D{v)^ 

(lib) X* = T + 5G{x, X — y)P{\v — w\){x — y) P iiD{x), 

(lie) = w — jR{y, y — x)P{\v — w\){w — v) P r]D{w), 

(lid) y* = y + SG{y, y - x)P{\v - w\){y - a;) + fiD{y). 

Herein, the functions P, P, and D play the same role as in the previous section, and are assumed 
to fulfil the assumptions introduced earlier. The constant parameters 7 G (0,1/2) and 6 G (0,1/2) 
control the ‘speed’ of attraction of opinions and repulsion of assertiveness, respectively. In addition, 
we introduce the function G which describes the in- and decrease of the individual assertiveness 
level. It is assumed to fulfil 

0 < G(t, X — y) < 1. 

Together, these assumptions guarantee that 1 ;*, re* G X and ^*, G X- The quantities 7 , y and /i, 
y are uncorrelated random variables with distribution 0 with variances and cr^, respectively, 
and zero mean, assuming values on a set S C M. 

To specify G, we make the assumption that highly assertive individuals gain more assertiveness 
and weakly assertive ones loose assertiveness in a collision. Therefore, we propose the following 
form: 

(12) G{x,x-y) = {l-x‘^)°^\x-y\, 

where the quadratic polynomial ensures that the assertiveness level remains in J. 


In this setting, we are led to study the evolution of the distribution function as a function 
depending on the assertiveness t G X, opinion w el and time t G M+, / = f{x^wR). Different 
to the previous section, we assume that the variation of the distribution f{x^wR) with respect 
to both variables, assertiveness x and opinion re, depends on collisions between individuals. The 
time-evolution of the distribution function / = f{x,wR) of individuals depending on assertiveness 
T G X, opinion w eX and time t G will obey a multi-dimensional homogeneous Boltzmann-type 
equation^ 


(13) 


^f{x,w,t) = P{f,f){x,W,t), 


where r is a suitable relaxation time which allows to control the interaction frequency. The 
Boltzmann-like collision operator Q which describes the change of density due to binary interac¬ 
tions is derived by standard methods of kinetic theory. The weak form of (13) is given by 


a Q{f, /) {x, w, t)4i{x, w) dw dx 

= ,W*) + (j){y*,v*) - (l>{x,w) - (p{y,v))f{x,v,t)f{y,w,t)dvdwdydx'^, 

for all test functions (j){x^w) where (•) denotes the operation of mean with respect to the random 
quantities 7 , 77 and /i, /i. 


3. Political segregation: ‘The Big Sort’ 

In this section we are interested in modelling the clustering of individuals who share similar 
political opinions, a process that has been observed in the USA over the last decades. 

In 2008, journalist Bill Bishop achieved popularity as the author of the book The Big Sort: 
Why the Clustering of Like-Minded Ameriea Is Tearing Us Apart [4]. Bishop’s thesis is that US 
citizens increasingly choose to live among politically like-minded neighbours. Based on county- 
level election results of US presidential elections in the past thirty years, he observed a doubling 
of so-called ‘landslide counties’, counties in which either candidate won or lost by 20 percentage 
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points or more. Such a segregation of political supporters may result in making political debates 
more bitter and hamper the political decision-making by consensus. 

Bishop’s findings were discussed and acclaimed in many newspapers and magazines, and former 
president Bill Clinton urged audiences to read the book. On the other hand his claims also met 
opposition from political sociologists [T] who argued that political segregation is only a by-product 
of the correlation of political opinions with other sociologic (cultural background, race, ...) and 
economic factors which drive citizens residential preferences. One could consider to include such 
additional factors in a generalised kinetic model, and potentially even couple the opinion dynamics 
with a kinetic model for wealth distribution mils]. 

‘The Big Sort’ lends itself as an ideal subject to study inhomogeneous kinetic models for opinion 
formation. In this context we are faced with spatial inhomogeneity rather than inhomogeneity in 
assertiveness as in Section In the following we propose a kinetic model for opinion formation 
when the individuals are driven towards others with a similar political opinion. 


3.1. An inhomogeneous Boltzmann-type equation. We study the evolution of the distri¬ 
bution function of political opinion as a function depending on three independent variables, the 
continuous political opinion variable w G [—1,1], the position t G on our (virtual) map and time 
t G / = f{x,w,t). In analogy with the classical kinetic theory of rarefied gases, we assume 
that the variation of the distribution /(x, re, t) with respect to the political opinion variable w de¬ 
pends on collisions between individuals, while the time change of distributions with respect to the 
position X depends on the transport term, which contains the ‘opinion-velocity field’ ^ = ^{x,w). 
The specific choice of $ is discussed further below. 

The exchange of opinion is modelled by binary collisions in the operator Q and has a similar 
structure as (11). Let x = {xi,X2) and y = (^1,^2) denote the positions of two individuals with 
opinion v and re, respectively. Then the interaction rule reads as: 


(14a) v"" = V — jK{\x — y\)P{\v — w\){v — re) + fjD{v), 

(14b) w"" = w — ^K{\y — x\)P{\v — re|)(re — r’) + r]D{w). 

In ( p^ we make the following assumptions: 

• Individuals exchange opinion if they are close to each other, i.e. their physical distance is 
less than a certain radius. This is modelled by the function K{z) = l^\z\<r} with some 
radius r > 0. 

• The functions P and D and the random quantities r), rj satisfy the same assumptions as 
in Section O 

We still need to specify the ‘opinion-velocity field’ 4>(T,re). As as first approach it helps to think 
of ^{x^ w) modelling the movement of individuals with respect to a given initial configuration: we 
consider and a number of counties Qi which constitute a disjoint cover of each with a county 
seat Ci = (c-^\cp^) G x [—1,1] with position in and an election result (by plurality voting 
system) which is either ‘red/Republican’ (re > 0) or ‘blue/Democratic’ (re < 0). In our numerical 
experiments presented later we will typically restrict ourselves to bounded domains C with 
no-flow boundary conditions. 

We assume that supporters of a party are attracted to counties which are controlled by the party 
they support. We model this effect by defining 4>(x,re) as a potential that drives the dynamics 
and is given by a superposition of (signed) Gaussians C{x) centered around with variance 
We assume also that stronger supporters, i.e. individuals with more extreme opinion values, are 
able to retain there positions. A possible way to take these effects into account is to choose 

(15) 4>(x,re) = sign(re)VC(T)(l — |t|^). 


The time-evolution of the distribution function / = f{x,w,t) of individuals with political 
opinion re G X at position x G at time t G will obey an inhomogeneous Boltzmann-type 
equation for the distribution function / = f{x^w,t) which is of the form 

+ diV:,(^{x,w)f) = ^Q(/,/)- 


(16) 
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In a two-party system as we consider it, quantities of interest are the supporters of the parties, 
given by the marginals 

(17) f]j(^x^t^ f(^x^w^t^ dw^ fji(^x^t^ i: /(t, re, t) dw. 

In simulations, elections could be run at time t by computing 

(18) Di{t) = f fD{x^t)dx= f f f{x,w,t)dwdx, 

JQi JQi J-1 

(19) Ri{t) = [ fR{x,t)dx= [ [ f{x,w,t)dwdx, 

JJ O j J 0 

(<2') _ 

to adapt the values of q ^ defining the attractive Gaussians C{x) in (15) accordingly in time. 


4. Fokker-Planck limits 

4.1. Fokker-Planck limit for the inhomogeneous Boltzmann equation. In general, equa¬ 
tions like ([T^ (and ( [Tq| )) are rather difficult to treat and it is usual in kinetic theory to study 
certain asymptotics, which frequently lead to simplified models of Fokker-Planck type. To this 
end, we study by formal asymptotics the quasi-invariant opinion limit ( 7 , ^ 0 while keeping 

cr ^/7 = A fixed), following the path laid out in [25] . 

Let us introduce some notation. First, consider test-functions 0 G C^’‘^([—1,1]) for some (5 > 0. 
We use the usual Holder norms 

|| 0 |U= ^ |T>||c + ^[i?>]co,., 

|a |<2 a =2 

where [h]co,5 = sup^^^ \h{v) — h{w)\/\v — w\^. Denoting by M.o{A), H C M the space of probabil¬ 
ity measures on H, we define 


Mp{A) = lee Mo 


f \r]\Pde{r]) 

JA 


< oc, p > 0 


the space of measures with finite pth momentum. In the following all our probability densities 
belong to M 2 -\-s and we assume that the density 0 is obtained from a random variable Y with 
zero mean and unit variance. We then obtain 


( 20 ) 


l^\vmv)dv = E[\ar,Y\P]=aP^E[\Yn 


where E[|F|^] is finite. The weak form of (16) is given by 

( 21 ) ^ f f{x^wA)<P{'^)dwdx-\- [ divx{^{x^w)f{x^wA))<P{'w)dwdx 

dt JxxR2 JXxR2 

= - [ Q{f,f){w)(j){w)dwdx 

T JIxR 2 


where 


I 


XxR2 


Q{fJ){w)(l){w)dw dx 


h \ / (0(re*)-h 0(r’*) — 0(r’) — (/)(x))/(x,x)/(x,R;) dxdredx) . 

2 \jx 2 xR 2 / 

To study the situation for large times, i.e. close to the steady state, we introduce for 7 <C 1 the 
transformation 

i='yt, x = 7a;, g{x, w, i) = f{x, w, t). 
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This implies f{x^w,0) = g{x^w,0) and the evolution of the scaled density g{x,w^i) follows (we 
immediately drop the tilde in the following and denote the rescaled variables simply by t and x) 

(22) ^ f g{x^w,t)(l){w) dw dxf div^{x^ w)g{x, w^ t))(l){w) dw dx 

= ^ / Q{g,g){w)(t){w)dwdx. 

1 '^ JXxR2 

Due to the collision rule Q, it holds 

— w = —7(7(x, p, w){w — x) + gD{w) 1. 


Taylor expansion of (j) up to second order around w of the right hand side of (22) leads to 

/— [ (p'(w) [—jC{x,y,v^w){w — v)gD{w)] g{x,w)g{x^v) dv dw dx\ 

JX2xR2 / 

+ /- [ (/)"(w) [—^C{x,y,v,w){w — v) d- gD{w)]‘^ g{x,w)g{x,v) dv dw dx\ 

\ 27 rJx 2 xR 2 / 

= — / (l)\w) [—jC{x,y,v,w){w — v))] g{x,w)g{x,v) dv dw dx 

7^ JX2xR2 

+ /I [yC{x^ y, p, w){w — x) + gD{w)] ‘^g{x^ w)g{x^ v) dv dw dx\ + cr^) 

\ 27 r 7x2 xR2 ' 

=- f (j)'{w)lC{x^w)g{x^w) dw dx 

^ 7xxR2 

H- [ ^‘^C‘^{x,y,v^w){w — v)‘^ d- jXD‘^{w) g{x,w)g{x^v) dv dw dx d- R{^^crr]), 

27 ^ 7x2 XR2 L J 

with w = KW* + (1 — i^)w for some k G [0,1] and 

= {d— j - (t>"{w)) 

\ 27 r 7x2 xR2 

X [— 7 ( 7 (x, y^ p, w){w — v) d- gD{w)]^ g{x^ w)g{x^ v) dv dw dx^ 

and 

IC{x,w) = J C{x,y,v,w){w — v)g{x,v) dv. 

Now we consider the limit 7 , cr^ ^ 0 while keeping A = cr ^/7 fixed. 

We first show that the remainder term i?( 7 ,cr^) vanishes is this limit, similar as in [25]. Note 
first that as 0 G 7 ^ 2 + 5 , by the collision rule @ and the definition of w we have 

\<p"{w) - riw)\ < U"\U\w - wf < U''\U\w* - w\^ = U"\\s \iC{x, y, v, w)iw -v) + yD{w)f . 
Thus we obtain 




L 


[—jC{x^y^v^w){w — v)d-gD{w)] g{x^w)g{x^v) dv dw dx). 
27 ^ WX2XR2 / 


Furthermore, we note that 

[gD{w) — 7(7(x, X, w){w — v)]^~^^ < {\jC{x^ y^ x, w){w — v)\‘^~^^ d- \gD{w)f~^^) 

^ 23 + 2 ( 5 | 7 | 2+(5 ^ 2 ^+^\rj\‘^+^, 

Here, we used the convexity of f{s) := and the fact that re,p G X and thus bounded. We 

conclude 

|i?(7, <^v)\ < ( 7 '+^ + ^ ( 7 '+^ ^ ^ X 
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Since 0 G and rj has variance we have (see ( 20 )) 


(23) 


\ri\^+^e{ri)dri 


= E 




2+5' 


= (A7)1+2e |y| 


■ I ‘2-\-5 


Thus, we conclude that the term i?( 7 , cr^) vanishes in the limit 7 , ^ 0 while keeping A = cr ^/7 

fixed. 


Then, in the same limit, the term on the right hand side of (22) converges to 
- f (l)\w)]C{x,w)g{x,w) dw dx ^f (l)'\w) [XD‘^{w)] g{x,w)g{x,v) dv dw dx 

^ JxxR2 2r Jx2xR2 

=- f (/)'{w)JC{x,w)g{x^w) dw dxf (j)" {w)M{x)D‘^{w)g{x^w) dw dx^ 

^ ixxR 2 2 r JxxR 2 

with M{x) = f g{x, v) dv being the mass of individuals in x. After integration by parts we obtain 
the right hand side of (the weak form of) the Fokker-Planck equation 

d 

(24) -^9{x,w,t) + diY^{^{x,w)g{x,w,t)) 

d X . XM(x) 52 

= ~ ) 4- 2^dw^ {D^{w)g{x,w,t)) , 


dw 


subject to no flux boundary conditions for the variable w (which result from the integration by 
parts). 

4.2. Fokker-Planck limit for the multidimensional model. We follow a similar approach as 
in the previous section, now in the two-dimensional setting of the opinion formation model ( 11 ) 
(cf. also [ 22 ]). Our aim is to study by formal asymptotics the quasi-invariant opinion limit where 
7 , J, ^ 0 while keeping cr ^/7 = cicr^/(c 2 ^) = A fixed, introducing the positive constants 

Cl = S/j and C 2 = cr^/cr^. 

We consider test-functions <p G C^’‘^([—1,1]^) for some (5 > 0. As above we use the usual Holder 
norms and consider denote by A4p(A), A C M the space of probability measures on A with finite 
pth momentum. In the following all our probability densities belong to At 2 +^ and we assume that 
the density 0 is obtained from a random variable Y with zero mean and unit variance. We then 
obtain 

|r?re(r?) dg = E[\aY\P] = <E[|Fn, WQifi) dfi = E[\aY\P] = af,E[\Yn 

where E[|T|^] is finite. The weak form of ( [l^ is given by 

^ f f{x^wA)4^{x^w)dwdx=— f Q{f ^ f){x^w)(j){x^w) dw dx 

JjxX XJjxX 


(25) 

where 


dt 


I ' 

Jjxl 


Q{f, /) {x, w)4>{x, w) dw dx 

f {(f>{x*,w*) + (f>{y*,v*)-(f>{x,w)-(p{y,v))f{x,v)f{x,w)dvdydwdx\. 

To study the situation for large times, i.e. close to the steady state, we introduce for 7 ^ 1 the 
transformation 

i = -ft, g{x, w, i) = f{x, w, t). 

This implies /(x,r;,0) = ^(x,r;, 0) and the evolution of the scaled density g{x,wA) follows (we 
immediately drop the tilde in the following and denote the rescaled time simply by t) 

(26) [ g{x^wA)4^{x^w)dwdx= — [ Q{g^g){x^w)(j){x^w) dw dx. 

JjxX IX JjxX 
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Due to the collision rules (11), it holds 

— w = —^R{x — y)P{\w — v\){w — v) r]D{w) <C 1, 

T* — T = SG{x — y)P{\v — w\){x — y) jiiD{x) <C 1. 

We now employ multidimensional Taylor expansion of 0 up to second order around (x^w) of the 
right hand side of (26). Recalling that the random quantities 77 , y have mean zero, variance 
and cr^, respectively, and are uncorrelated, we can follow along the lines of the computations of 
the previous subsection, to obtain 

— [ ^^{x^w)[—jR{x — y)P{\w — v\){w — v)-\-r]D{w)]g{x^w)g{x^v)dvdydwdx 

fT 


+ 


f ^ 

7^ Jj 


Jj^xx^ dx 
d‘^(j) 


(x, w) 6G{x — y)P{\v — w\){x — y) P yD{x) 


X g{x, w)g{x^ v) dv dy dw dx 
+ (-7 f Pt(^x,w)[-^R{x-y)P{\w-v\){w-v)+r]D{w)f 

X g{x^ w)g{x, v) dv dy dw dx 


P 


/J_ / ^ 

T'T Jn 


+ 


\ 2^t 

tJ. 


J^xX‘^ 


dx‘^ 


(x, w) 5G{x — y)P{\v — w\){x — y) P yD{x) 


1 2 


X g{x^ w)g{x^ v) dv dy dw dx 


5^0 
Jj 2 y^x‘^ dwdx 


(x, w) [—^R{x — y)P{\w — v\){w — u) + gD{w)] 


X 6 G{x — y)P{\v — w\){x — y) P yD{x) g{x^w)g{x^v) dv dy dw dx 

+ -R(7,crr?,0'^) 

1 f dcj) S f dcj) 

=-/ ^—(x^w)JC(x^w)g(x^w)dydwdx -/ -^ix^w)Cix^w)gix^w) dv dw dx 

X Jxxj^ dw Jx2^j dx 


<J^ 


hL 


2it Jx2^j 2 dw^ 


(x, w) — y)P^{\w — v\){w — vy P a^D^{w) 


X g{x,w)g{x,v) dv dy dw dx 


-[ 

Jx2 


d‘^<p 


+ -— / --^{x,w) 6 G {x-y)P {\v-w\){x-y) P cr D (x) 


P 


27 ^ Jx^xj^ dx^ 
d^(j) 


X g{x^w)g{x^v) dv dy dw dx 


-[ 

Jx 


2^t Jx‘^xj^ dwdx 


(t, w) ^SR{x — y)P^{\w — v\){w — v)G{x — y){x — y) 


X g{x,w)g{x^v) dv dy dw dx 
+ R( 7 ,cr^,cr^), 


where 

(27a) IC{x,w) = j R{x — y)P{\w — v\){w — v)g{x^v) dv^ 

(27b) C{x,w)= [ G{x-y)P{\v-w\){x-y)g{x,v)dy. 

Jj 

and R( 7 , cr^, cr^) is the remainder term. Our aim is now to consider the formal limit 7 , (5, cr^, cr^ ^ 0 
while keeping cr ^/7 = cicr^/(c 2 ^) = A fixed, recalling ci = 5/^ and C 2 = cf^IcFt)- The remainder 
term ^( 7 , cr^,cr^) depends on the higher moments of the (uncorrelated) random quantities and 
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can be shown to vanish in this limit, similar as in the previous subsection (we omit the details). 
In the same limit, the term on the right hand side of (26) then converges to 


^ Jx 


XxJ‘^ 


(x^w)]C(x^w)g(x^w) dy dw dx — — [ ^(x^w)C(x^w)g(x^w)dvdwdx 

^ Jx^xj dx 

hi 

- / -^{x,w)IC{x,w)g{x^w) dy dw dx -/ -^{x,w)C{x,w)g{x,w) dv dw dx 

^ Jxxj^ ^ Jx'^xj 

L 


^^xj 

(j) 

w)XD^(w)g(x^ w)gix^ v) dv dy dw dx 

Jx^xj^ 

0^ (j) 

—-^(x, w)c2XD‘^{x)g{x, w)g{x, v) dv dy dw dx 
X2xj^2 dx 


X 


d^(\) 

}XxJ‘^ 


(t, w)M{x)D‘^ (w)g{x^ w) dy dw dx 
d^(j) 


+ ^ [ ^^{x,w)M{x)D‘^{x)g{x,w)dydwdx, 
2r Jxxj^ dx^ 


with M{x) = f g{x, v) dv being the mass of individuals in x. After integration by parts we obtain 
the right hand side of (the weak form of) the Fokker-Planck equation 


(28) 


' ^ {D^{w)g{x,w,t)) + {D^{x)g{x,w,t)) , 


2r dw‘^ ^ / 2r dx^ 

subject to no flux boundary conditions which result from the integration by parts, with the nonlocal 
operators /C, jC given in (27). 


5. Numerical experiments 


In this section, we illustrate the behaviour of the different kinetic models and the limiting 
Fokker-Planck-type equations with various simulations. We first discuss the numerical discretisa¬ 
tion of the different Boltzmann-type equations and the corresponding Fokker-Planck-type equa¬ 
tions, and then present results of numerical experiments. 


While the multi-dimensional Boltzmann-type equation (13) can be solved by a classical kinetic 


Monte Carlo method, the Monte Carlo simulations for the inhomogeneous Boltzmann-type equa¬ 
tions ^ are more involved. On the macroscopic level the high-dimensionality poses a significant 
challenge - hence we propose a time splitting as well as a finite element discretisation with mass 
lumping in space. 


5.1. Monte Carlo simulations for the multi-dimensional Boltzmann equation. We per¬ 
form a series of kinetic Monte Carlo simulations for the Boltzmann-type models presented in 
Sections [H and |2|^ In this kind of simulation, known as direct simulation Monte Carlo 
(DSMC) or Bird’s scheme, pairs of individuals are randomly and non-exclusively selected for 
binary collisions, and exchange opinion (and assertiveness in the multi-dimensional model from 
Section 2|2.2 ), according to the relevant interaction rule. 

In each simulation we consider N = 5000 individuals, which are uniformly distributed in X x J" 
at time t = 0. One time step in our simulation corresponds to N interactions. The average steady 
state opinion distribution / = f{x,w,t) is calculated using M = 10 realisations. To compute a 
good approximation of the steady state, each realisation is carried out for n = 2 x 10^ time steps, 
then the particle distribution is averaged over 5000 time steps. The random variables are chosen 
such that gi and gi assumes only values u = ±0.02 with equal probability. We assume that the 
diffusion has the form D{w) = (1 — w‘^)^ to ensure that the opinion w remains inside the interval 
X. The parameter a is set to g = 2 if not mentioned otherwise. 
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5.2. Probabilistic Monte Carlo simulations for the inhomogeneous Boltzmann equa¬ 
tion. The Monte Carlo simulations for the inhomogeneous Boltzmann-type equations § are more 
involved. In equation (§ the advection term is of conservative form, hence the transportation step 
cannot be translated directly to the particle simulation. Pareschi and Seaid [21] as well as Herty et 
al [IHI propose a Monte Carlo method, that is based on a relaxation approximation of conservation 
laws. It corresponds to a semilinear system with linear characteristic variables. We shall briefly 
review the underlying idea for the conservative transportation operator in (§ in the following. 

Let us consider the linear conservation law in one spatial dimension for the function p = p{x^ t) 
with a given flux function 6(x, t) : M x [0, T] ^ M: 


(29) 


dtP + d:^(b{x, t)p) = 0, p{x, 0) = po{x). 


Note that ( |29| ) corresponds to the transportation step in a splitting scheme in a Monte Carlo 
simulation for Equation ( [2^ can be approximated by the following semilinear relaxation 

system: 


(30) 


dtp + dxV = 0 , dtv + bdxU ■■ 


-i(n. 


b{x,t)p) 


with initial conditions p(t, 0) = po{x) and p(x,0) = b{x,0)po{x) and b^e > 0. The function v 
approaches the solution at local equilibrium v = (p{u), if —b < b{x,t) < b for all x. Note that 
system (30) has two characteristic variables given by p ± Vb, which correspond to particles either 
moving to the left or the right with speed Vb. Hence, we introduce the kinetic variables p and q 
with 


p = p ^ q and v = b{p — q). 

Then the relaxation system can be written as follows: 

The numerical solver is based on a splitting algorithm for system It corresponds to first 

solving the transportation problem and then the relaxation step. While the transportation step is 
straight forward, the relaxation step can be interpreted as the evolution of a probability density. 
Since 

P>0, ^^>0, --h- = l and dtp = 0 

P P 

in the relaxation step we can calculate the solution explicitly and obtain 

\ t 1,, li 1., 

p =-e^ -\ — =b(x,t)p and q = —- =b(x,t). 

^2 2b ^ ^ 2 2b ^ ^ 

Then the solution at the time = (n + l)At reads as: 

p(x, = (1 — A)p(t, f^) + A6(x, f^)p{x^ g(T, = (1 — A)g(T, f^) — A6(x, f^)p{x^ 

with A = 1 — Let p'^ = p(t, q^ = g(x, and p'^ = p^ ^ q^. Since 0 < A < 1 we can 

define the probability density 

{ pn jpU if ^ = y/b^ 

q^/p^ if^ = -A 

0 elsewhere. 

Since 0 < < 1 and relaxation step can be interpreted as the evolution of 

a probability function 

pn+l(^)^(l_A)pn(^)^AP-(0, 

where E^{^) = b{x,E)p'^ if ^ and = —b{x,E)p'^ if ^ = —V^. So in the probabilistic 

Monte Carlo method a particle either moves to the right or the left with maximum speed \/b in 
the transportation step. Then the two groups are resampled according to their distribution with 
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respect to the equilibrium solution in the relaxation step. For further details we refer 

to [H]. 


5.3. Fokker-Planck simulations. The discretisation of the Fokker-Planck-type equation (24) is 
based on a time splitting algorithm. The splitting strategy allows us to consider the interactions 
in the opinion variable, i.e. the right hand side of equation (24), and the transport step in space 
separately. 

Let At denote the size of each time step and = kAt. Then the splitting scheme consists of a 
transport step S^{g, At) for a small time interval At: 

— {x,w,t) + div^ {(f>{x,w)g*{x,w,t)) = 0 , 


g*{x,w,0) =go{x,w), 


(32a) 

(32b) 

and an interaction step S^{g, At): 

tlC{x,w,t)g^{x,w,t)) + ^^{D^{w)g''{x,w,t)), 

g^{x,w,0) = g*{x,w,At). 


(33.) 

(33b) 

The approximate solution at time t = is given by 

g'^+\x,w) = S 2 {g*’^+\At/ 2 ) o At) o 52 (/, At/ 2 ), 

where the superscript indizes denote the solution of and at the discrete time steps t^ = kAt 
and t ^+2 = (/c + ^)At. Both equations are solved using an explicit Euler scheme in time and a 
conforming finite element discretisation with linear basis functions i = 1 ,... (#points), in space 
as well as opinion. Then the discrete transportation step ( [3^ in space has the form 

M 


■) - gj^ix, •) = At(T[g!^{x, •)]) 


where M = {pi^Pj)ij corresponds to the mass matrix for element-wise linear basis functions and 
T = {^{x)\/pi^pj)ij denotes the matrix corresponding to the convective field 4> of the form 
In the interaction step we approximate the mass matrix by the corresponding lumped mass matrix 
M, which can be inverted explicitly. The discrete interaction step (|3^ reads as: 

M 


'gp\;w) - 9 H;U>)] = ^t(tc[K[gf],gf{;w)] + ^^L[gf{;w)]), 


with the discrete convolution-transportation matrix C = {K{x^w)pi^Vpj) and Laplacian L = 
{D{w)Vpi^Vpj). Here the vector K corresponds to the discrete convolution operator /C, which 
has been approximated by the midpoint rule: 

K{x,wi) = '^C{x,Vk,wi){wi - Vk)g{x,Vk)Aw. 
k 

5.4. Numerical results: opinion formation and opinion leadership. We present numerical 
results for the kinetic models for opinion formation including the assertiveness of individuals from 
SectionWe consider the inhomogeneous Boltzmann-type model of Section 2 | 2.1 and the multi¬ 
dimensional Boltzmann-type equation of Section [ 2 |^ which are discretised using the Monte Carlo 
methods presented in the previous sections. 


5.4.1. Example 1: Influence of the interaction radius r. First we would like to illustrate the in¬ 
fluence of the interaction radius r. We assume that the functions in both models which relate to 
the increase or decrease of the assertiveness (see and (El) define a constant increase of the 
individual assertiveness level, i.e. 

G(t, w) = 1 and G(t, x — y) = {1 — x‘^)^ tanh (^k{x — y )), 

with k = 10. We start each Monte Carlo simulation with an equally distributed number of 
individuals within (u,t) G (—1,1) x (— 1 , 1 ). We expect the formation of one or several peaks 
at the highest assertiveness level due to the constant increase caused by the functions G and G. 
Figure which shows the averaged steady state densities for = 7 = 0.25 and two different radii. 
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r = 0.5 and r = 1. For r = 1 we observe the formation of a single peak located at the highest 
assertiveness level, re = 1, in Figure in the inhomogeneous and multi-dimensional model. In the 
case of the smaller interaction radius, r = 0.5, two peaks form at the highest assertiveness level. 
These results are in accordance with numerical simulations of the original Toscani model 


0.6 J-' 0.35- 



^ ^ opinion ^ ^ opinion 

assertiveness assertiveness 

(a) Inhomog. r = 1 (b) Inhomog. r = 0.5 



(c) Multi-dim. r = 1 


(d) Multi-dim r = 0.5 


Figure 2. Example 1: Steady state densities of the inhomogeneous and multidi¬ 
mensional Boltzmann-type opinion formation model for different interaction radii 

r. 


5.4.2. Example 2: Different ehoiees of G and G. Next we focus on the influence of the functions 
G and G, which model the increase of the assertiveness level in either model (see § and in)- 
We choose the same initial distribution as in Example 1 and set 

G(t, w) = tanh(A:T) and G(t, x — y) = {1 — x‘^)^\x — y\. 

Both functions are based on the following assumption: individuals with a high assertiveness level 
gain confidence, while those with a lower level loose. We expect the formation of peaks close to the 
highest and lowest assertiveness level. This assumption is confirmed by our numerical experiments, 
see Eigurej^ We observe the formation of peaks close to the maximum and minimum assertiveness 
level; the number of peaks again depends on the interaction radius r. 
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(c) Multi-dim. r = 1 


(d) Multi-dim r = 0.5 


Figure 3. Example 2: Steady state densities of the inhomogeneous and mul¬ 
tidimensional Boltzmann-type opinion formation model for different interaction 
radii. 


5.4.3. Example 3: Preferred opinions promoting leadership. Our last example illustrates the rich 
behaviour of the proposed models. We assume that leadership qualities are directly related to the 
opinion of an individual. Individuals with a popular, ’mainstream’ opinion, i.e. w = ±0.5 gain 
confidence, while extreme opinions are not promoting leadership. To model this we set 

G{x,w)=Gix,w,x-y) = —± 

V27rcr V / 

with C7^ = 1/8 and 


(34) 


i?(a 


- 2/) = (^ - ^ tanh {k{x - y 




with k = 10. The second factor in (34) models the assumption that individuals change their 
opinion more if they have a low assertiveness level. 

At time t = 0 we equally distribute the individuals within (re, x) G [—0.25, 0.25] x [—0.75, —0.25] 
- hence we assume that initially no extreme opinions exist and no leaders are present. The 
interaction radius is set to r = 1. Figure [^illustrates the very interesting behaviour in this case. 
In the multidimensional simulation we observe the formation of a single peak at a low assertiveness 
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level, while individuals with a higher assertiveness level group around the ‘mainstream’ opinion 
w = ±0.5. In the case of the inhomogeneous model, the potential ^ initiates an increase of the 
assertiveness level for all individuals. Therefore we do not observe the formation of a single peak 
at a low assertiveness level in this case. 



(a) Inhomog. (b) Multi-dim. 

Figure 4. Example 3: Steady state densities of the inhomogeneous and multi¬ 
dimensional Boltzmann-type opinion formation model. 


5.5. Numerical results: ‘The Big Sort’ in Arizona. In our final example we present numer¬ 
ical simulations of the corresponding Fokker-Planck equation (24), which illustrate ‘The Big Sort’ 
by considering the state of Arizona. Arizona is a state in the southwest of the United States with 
15 electoral counties. In recent years the Republican Party dominated Arizonas politics, see for 
example the outcome of the presidential elections from the years 1992 to 2004 in Figure The 
colours red and blue correspond to Republicans and Democrats, respectively. The colour intensity 
reflects the election outcome in percent, i.e. dark blue corresponds to Democrats 60-70%, medium 
blue to Democrats 50-60% and light blue to Democrats 40-50%. Similar colour codes are used 
for the Republicans. The election results illustrate the clustering trend as the election results per 
county become more and more pronounced over the years. 

We solve the Fokker-Planck equation (24) on a bounded physical domain U C corresponding 
to the state Arizona, divided into 15 electoral counties, see Figure |6 (a) | We employ the numerical 
strategy outlined in Section 5|5.3 We discretise the physical domain in 9356 triangles, the opinion 
domain J' = [—1,1] into 100 elements. The time steps are set to At = 2 x 10“^. 

We choose an initial distribution which is proportional to the election result in 1992, 


/(t,r;,0) = 0.5 ± 


w{l - w)fD,1992{x) 

w{l^w)fR^l992{x) 


for re > 0, 
for re < 0, 


where fD,i 992 and fn ,1992 correspond to the distribution of Democrats and Republicans estimated 
from the election results in 1992. We approximate the distributions by assigning different constants 
to the respective percentages in the electoral vote, i.e. 


(35) 


fD,R{x) = < 


0.25 

0.5 


0.75 


if the election outcome is < 40 — 50% 
if the election outcome is between 50 — 60% 
if the election result is > 60%. 


The potential 4>, given by (15), attracts individuals to the counties controlled by the party they 
support. We assume that the potential C = C(x) is directly related to the electoral results of the 
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(a) 1992 


(b) 1996 


(c) 2000 



(d) 2004 


Figure 5. Results of the presidential elections in Arizona [19]. The colour in¬ 
tensities reflect the election outcome in percent, i.e. dark blue (red) to Democrats 
(Republicans) 60-70%, medium blue (red) to Democrats (Republicans) 50-60% 
and light blue (red) to Democrats (Republicans) 40-50%. 



(a) Computational domain 
Q 



Figure 6. Computational domain O representing the 15 electoral counties of 
Arizona and potential C corresponding to the election results in 1996. 


year 1996, and satisfies the following PDF: 

C{x) + sAC{x) = fi 996 {x) for all x G 11, VC(t) • n = 0 for all x G 90. 


The Laplacian is added to smooth the potential C, the right hand side corresponds to the election 
results in the year 1996 (using the same constants as in (35) with ± sign corresponding to the 
Republicans and Democrats respectively). Note that we choose Neumann boundary conditions to 
ensure that individuals stay inside the physical domain. The calculated potential C is depicted in 
Figure [6(1^ The simulation parameters are set to A = 0.1, r = 0.25 and the interaction radius to 
r = 5. 


Figure shows the spatial distribution of the marginals ([ 13 , i.e. 

3 d 

f{x,w,t)dw and fR{x,t) = / f{x,w,t)dx, 

1 Jo 

which correspond to Democrats and Republican, respectively, at time t = 0.5. We observe the 
formation of two larger clusters of democratic voters in the south and the Northeast of Arizona. 
The republicans move towards the Northwest as well as the Southeast. This simulation results 
reproduce the patterns of the electoral results from 2004 fairly well, except for the second county 
from the right in the Northeast (Navajo). However, given the electoral data from 1992 and 
1996 only it is not possible to initiate such opinion dynamics. This would require more in-depth 


fD{x,t) 


r 
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(a) Distribution of democrats 


(b) Distribution of republicans 


Figure 7. The Big Sort: distribution of Democrats, foix, t) = /(x, w, t) dw, 

and Republicans, fR{x, t) = f(x,w, t) dw. 


information such as the demographic distribution within the counties or the incorporation of 
several sets of electoral results. We leave such extensions for future research. 

6. Conclusion 

We proposed and examined different inhomogeneous kinetic models for opinion formation, when 
the opinion formation process depends on an additional independent variable. Examples included 
opinion dynamics under the effect of opinion leadership and opinion dynamics modelling political 
segregation. Starting from microscopic opinion consensus dynamics we derived Boltzmann-type 
equations for the opinion distribution. In a quasi-invariant opinion limit they can be approximated 
by macroscopic Fokker-Planck-type equations. We presented numerical experiments to illustrate 
the models’ rich behaviour. Using presidential election results in the state of Arizona, we showed 
an example modelling the process of political segregation in the ‘The Big Sort’, the process of 
clustering of individuals who share similar political opinions. 
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